arXiv:1505.04109vl [cond-mat.stat-mech] 15 May 2015 


Condensation in models with factorized and 
pair-factorized steady states 


M. R. Evans and B. Waclaw 

SUPA, School of Physics and Astronomy, University of Edinburgh, Peter Guthrie 
Tait Road, Edinburgh EH9 3FD, United Kingdom 

E-mail: martinSph. ed. ac. uk, bwaclawOstaf fmail .ed.ac.uk 

Abstract. Non-equilibrium real-space condensation is a phenomenon in which a 
finite fraction of some conserved quantity (mass, particles, etc.) becomes spatially 
localised. We review two popular stochastic models of hopping particles that lead to 
condensation and whose steady states assume a factorized form: the zero-range process 
and the misanthrope process, and their various modifications. We also introduce a new 
model - a misanthrope process with parallel dynamics - that exhibits condensation and 
has a pair-factorized steady state. 


PACS numbers: 89.75.Fb, 05.40.-a, 64.60.Ak 

Date May 18, 2015 


Condensation in models with faetorized and pair-factorized steady states 

1 . Introduction 


2 


Real space condensation is a non equilibrium phase transition which occurs in various 
contexts such as granular clustering, traffic jams, wealth condensation or simulations of 
polydisperse hard spheres [H El [31 S]. In all of these systems there is some conserved 
quantity (mass, wealth or volume for example) which is transported through the system. 
If the global density of this quantity is above a critical value, a hnite fraction condenses 
onto a single lattice site or localized region of space. 

Surprisingly, many features of condensation are captured within a simple lattice 
model known as the zero-range process (ZRP) (for review see |5l |6] ). In the simplest 
one-dimensional asymmetric version of this model, a particle moves from site i to i-fl of 
a one-dimensional periodic lattice with rate u{m), where m is the occupancy (number 
of particles) of the departure site i. As the rates are totally asymmetric a current 
always flows and detailed balance cannot be satished, thus the stationary state is non- 
equilibrium. The great advantage of this model is that its non-equilibrium stationary 
state has a simple factorized form which is amenable to exact analysis. The structure 
for P{{mi}), the probability that each site 1 < i < L contains mass is 

1 ^ 

^ n rn^ ■ (1) 

^ i=l 

Thus the numerator in ([T]) contains one (non-negative) factor /(m*) for each site i and 
f{m) is known as the single-site weight function and depends on u{m) as 

m ^ 

f(m) = IT —— for m > 1 and /(O) = 1. (2) 

1 u{i) 

2=1 ^ ^ 

The denominator Zl is the normalization or nonequilibrium partition function 

oo L 

Zl= (3) 

{ mi = 0 } *=1 


In ([T]) and ([3]) the Kronecker delta imposes the constraint that the total mass (or number 
of particles) M in the system is conserved. 

It turns out that the zero-range process is not the only model leading to the 
factorized steady state ([T]). In recent years, it has been shown that other models in 
which the hop rate depends not only on the occupation of the departure site but also on 
other variables can also exhibit steady states that factorize exactly [aiHiEiiiniiiiiiniiis] 
or approximately [H] over the sites of the system. It has also been discovered that there 
exist various classes of models in which the steady state factorizes over pairs of sites 

[laiieiin]. 




( 4 ) 


2 = 1 


In dl]) 5 f(mj,mj+i) is the pairwise weight. In the following we will hrst briefly review 
some of the models with fully factorized and pair-factorized steady states and discuss 
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conditions under which condensation occurs in these models. We shall then introduce 
a new model, a discrete time variant of the misanthrope process [18], with a pair- 
factorized steady state in which the hopping rate depends on the occupations of the 
departure and arrival sites. Such discrete time schemes are often used in the simulation 
of traffic ffow and pedestrian dynamics. In the zero-range process it has been shown 
that a generalisation to discrete time dynamics still results in a factorized stationary 
state [19]. In this work we show that the misanthrope process may have a pair-factorized 
stationary state under discrete time dynamics and we establish conditions under which 
this holds. Further, we present an analysis of condensation in such pair factorized states. 


2. Zero-range process 


The zero-range process is specihed by the hop rate u{m) which determines the properties 
of the steady state through the single-site weight f{m) from Eq. (I2|). It is important 
to note that any exponential factor Aq"^ in f{m) does not change the steady-state 
properties since it appears in Eq. ([T]) as a constant prefactor = A^q^ due to 

the hxed total mass M and number of sites L. From now on we will generally suppress 
such exponential factors in f{m). 

Condensation in models with factorized stationary states occurs in the limit L —)■ oo 
and hxed density of particles p = M/L when the asymptotic (large m) behaviour of f{m) 
(modulo any exponential factors) is the following: 

I f{m) m with 7 > 2. The critical mass density above which condensation occurs 

is hnite but its numerical value depends on the particular form of f{m) and not 
only on its asymptotic behaviour. The fraction p/pc — 1 of all particles goes into 
the condensate. We will refer to this behaviour as standard condensation. 0 

II f{m) increases with m more quickly than exponentially, e.g., as ~ m!. This leads 

to so called strong (or complete) condensation - the critical density pc = 0 and a 
fraction of particles tending to one in the thermodynamic limit is located at one 
site. 


It can be shown that standard condensation (I) occurs in the ZRP when the hop rates 
in the limit of large m asymptotically approach some positive value ft as 


u{m) 


1 + 


7 


m 


(5) 


with 7 > 2 , or more slowly than 1/m. On the other hand strong condensation occurs 
when u{m) —?■ 0 as m ^ oo. For example, u{m) = 1/m yields /(m) = m!. 

To see why the condensation happens in the two generic cases highlighted above, 
we shall follow a standard approach [5]. Treating the steady-state probability as the 


f Actually condensation also occurs if /(m) decays more quickly than l/m'*', e.g., as a stretched 
exponential. 
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statistical weight of a given configuration, and defining the grand-canonical partition 
function 

L 

G(^) = E z ^', mi) = En f{m,)z^^ = F{z)\ (6) 

{mi} {rrii} i=l 

where 

OO 

F(z) = E/(™)^”. (7) 

m=0 

we see that the phase transition, signaled by a singularity of G{z) at some Zc, is possible 
only if the series F{z) has a finite radius of convergence Moreover, the density 
calculated as a function of fugacity z from the grand-canonical partition function; 

F'(z) 

( 8 ) 

must yield a value Pc < oo as z Zc so that the singularity in G{z) and accompanying 
phase transition occur at finite density. This is only possible if either f{m) decays as a 
power law in which case we may have a finite pc (case I) or f{m) grows very fast with 
m in which case Zc = 0 and Pc = 0 (case II). 

Thus the grand canonical ensemble can only realise densities p < p^. When p > Pc 
one must work in the canonical ensemble (fixed number of particles) [51 EO]- If turns 
out that the excess mass Mgx = M — Lp^ condenses onto a randomly selected lattice 
site and forms the condensate. The remainder of the system (referred to as the fluid) is 
described by the grand canonical ensemble at the critical density pc [211 ESI EQl ESI EH. 

3. Generalized Class of Models with Factorized Stationary State 

So far we have discussed the ZRP as an example of a simple model with a factorized 
stationary state o. More generally one can ask, when does a stochastic mass transport 
model have a such a factorized stationary state if the hop rate depends only on the 
state of the departure site? To this end a class of models was studied in [19] that 
generalises the ZRP in a number of different ways while maintaining the factorization 
of the steady state. First, more than one unit of mass can be transferred from site i to 
i + 1. Second the dynamics consists of discrete time update and simultaneous transport 
of mass at different locations is possible at an update. More precisely, in each time 
step, some number 0 < /ij < m* of particles depart from site i and move to site i + 1 
with probability 0(/ij|mj) which is known as the chipping kernel. For conservation of 
probability we require = 1- 

It was shown [19] that a necessary and sufficient condition for the stationary state 
of this class of models to factorize is that the chipping kernel takes the form 

u{p)v{m — p) 
f{m) 


0(/i|m) 


(9) 
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Figure 1. Definition of the Misanthrope process: a particle hops from site with m to 
site with n particles with rate u(m,n). 


where the single site weight f{m) is given by 

m 

f{m) = ^u{iJ,)v{m-la) . (10) 

1 - 1=0 

In these expressions u and v are positive functions. Expression ([9]) implies a factorization 
of the chipping kernel into a factor u{fi) which depends on the mass transferred and a 
factor v{m — fi) which depends on the mass which remains. The model can be further 
generalized to a continuous mass variable |19| but we shall not consider that here. The 
model was also considered on an arbitrary graph rather than a periodiic chain and a 
condition similar to (j9]) was show to be sufficient for factorization |25| . 


4. Misanthrope process 


As already stressed the dehning feature of the ZRP and the class of models just discussed 
is that transition rates or probabilities for the transfer of mass between sites depend only 
on the departure site and not on the destination site. It is natural to consider more 
general models in which for example a hop rate takes the form u{m,n) where m is the 
occupancy of the departure site i and n is the occupancy of the destination site z + 1. 
This type of model is variously referred to as a misanthrope or migration process. 

We now dehne the model that we consider. As in the ZRP case, M particles reside 
on sites of a ID closed chain of length L; each site i carries m* particles, and the 
conservation of particles requires that ~ The only difference with the ZRP 

is that a particle hops from site i to site z + 1 with rate zz(mj, mj+i) which depends on 
the occupancies of both the departure and the arrival site, see Fig. [1] 

This model has a factorized stationary state ([I]) when certain conditions on u{m, n) 
are satished. Here, we simply quote the constraint on u{m,n) without proof: 


u{m, n) = u{m + 1, n — 1)- 


zz(l, m)u{n, 0) 


u{m + 1, 0)zz(l, rz — 1) 

It can be shown that this relation actually reduces to two conditions; 


+ u{m, 0) — u{n, 0). (11) 


u{n, m) 
u{n, m) 


. ^ zz(l, m)zz(rz, 0) 

uim + 1, rz — 1)—--—;- 

u{m + 1, 0)zz(l, rz — 1) 

zz(m, rz) = zz(rz, 0) — zz(rrz, 0). 


( 12 ) 

(13) 
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These conditions were first written down by Cocozza-Thivent [TH]. Under these 
conditions the single site weights f{m) obey a recursion 

f(n) = fin - (14) 

Conditions fll2lll3p uniquely define all rates u{n,m) and the single-site weights f{m) in 
terms of a set of basic hop rates which we denote 


u{m,0) = ym, (15) 

u{l,n) = Xn. (16) 


Iterating Eq. flTT)) and using the definitions of ym, Xn fll5lll6p we obtain the following 
expression for the single-site weight /(n): 


fin) 



(17) 


in which, as noted previously, we may suppress the exponential factor /(O) • 

Equations flT^ and m can be rewritten as two recursion relations which allow 
one to find u{m, n) for m > 2, n > 1: 

u{m -I- 1, n — 1) = u{n, , (18) 

Vn 

u{n,m) = u{m,n)-ym + yn- (19) 


By iterating these equations one obtains unique expressions for all u{m,n). However, 
there is an additional condition that u{m, n) should be non-negative for all m, n for it to 
be a hopping rate. This imposes some constraints on ym,Xn which cannot be expressed 
in a closed form. Therefore it remains an open problem to determine precisely which Xn 
and y^n lead to a physical model with non negative hopping rates. However, there exists 
a special case of u{m, n) for which more progress has been made and we shall discuss it 
now. 


5. Factorized hop rate in the misanthrope process 

In recent work [26l |27] we considered the special case 

u{m,n) = w{m)v{n) (20) 

which corresponds to hopping rates whose dependence on departure and destination site 
factorizes. One can check that for this form of the hop rate equation flT^ is automatically 
fulhlled. When v{n) = const, equation ffT3|) is also fulfilled and we recover the ZRP with 
u{ni,n) = w{m). When v{n) ^ const, equation (1T3)) leads to the relation between w{n) 
and v{n): 

tc(n) = U[n(n) — n(0)], (21) 

with some arbitrary, non-zero constant C. Thus 
u{m,n) = C[v{m) — n(0)]n(n) 


(22) 
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is the form of a factorized hop rate that yields a factorized steady state ([T]) with the 
single-site weight given by 


/(m) = Yl 


v[i 


1 ) 


2=1 


v{i) — n(0) 


(23) 


Before we discuss the condition for condensation in this model, we shall briefly review 
two simple choices of v{m) that lead to previously studied models. 


5.1. Partial Exclusion 

Our first simple example of the factorized hopping rate (12^ is O = — 1 and v{n) = N — n 
with some integer iV > 0 we have 

u{m, n) = m{N — n) . (24) 

If iV = 1 this reduces to the asymmetric simple exclusion process where the occupancy 
of each site is limited to 1 and m( 1, 1) = 0. Similarly, the case of general integer 
TV > 0 corresponds to ‘partial exclusion’ |28| where each site of a lattice contains 
at most N particles. In this context the rate fl2H) may be understood as each of m 
particles attempting hops forward to the next site with rate one and the hopping attempt 
succeeding with probability N — n where n is the occupancy of the destination site. 


5.2. Inclusion Process 
If we take C = 1 and 

v{n) = n + d (25) 

in fl22|) where d is a positive constant, we obtain 

u{m,n) = m{n + d) . (26) 

This is the hop rate for the so-called Inclusion Process studied in [21]. In the limit d —)■ 0 
this model exhibits a distinct form of condensation. 


6. Condensation in the misanthrope process with factorized hop rates 

We are interested in a factorized form of u{m, n) from Eq. (j22|) such that Eq. (lT4j) gives 
f{n) n with 7 > 2 which as we know leads to standard condensation. It turns out 
that the standard condensation can occur through two contrasting types of dynamics. 
One mechanism is through the hops rates decaying sufficiently slowly with n, m. This 
can be achieved through 

w(m) =/T ("l - —) (27) 

V mJ 

which leads to 

o/T(n(0) - jd) 


u{m, n) = /5(n(0) — (3) 


n 


m 


(28) 
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which decays with m in a similar fashion to the ZRP case. The other mechanism is for 
the hop rates to increase with m, and n as 

u{m,n) ~ {mny (29) 

with 7 > 2, which is equivalent to 

v{m) ~ m'^. (30) 


We refer to the latter case as explosive condensation. Explosive condensation exhibits 
strikingly different dynamical properties to the ZRP-like condensation. In particular, 
the condensate emerges on a time scale which vanishes with system size L as (InL)i-^ 
for 7 > 2 |26| . This is in contrast to the case (|27|) or the ZRP case for which the time 
increases with L as ~ L^. 

Interestingly, for the misanthrope process the existence of condensation depends not 
only on the asymptotic behaviour of v{m) but also on n(0). This should be contrasted 
with the ZRP for which it is only the asymptotic decay of the hop rate that determines 
condensation. To illustrate this point consider the case EH 


n(0) < 1, n(m) = 1 + 


m + 1 


(31) 


In this case one can find closed form expressions for the weights f{n) and generating 
function F{z) and one may determine the critical density given hy z ^ 1: 

^ 3n(0)-2 ^ ^ 

We see that Pc 0 as n(0) —?■ 1 and pc —t cx) as n(0) —)■ 2/3. Conseequently 
for n(0) < 2/3 there is no condensation, but for 2/3 < n(0) < 1 there is standard 
condensation and for n(0) = 1 there is strong condensation, even though v{m) is the 
same in all cases for m > 0. 


7. Pair-factorized steady states 


So far we have discussed the processes in which the stationary probability P{mi ,..., mi) 
factorizes over sites of a Id closed chain. One can consider generalisations of this 
structure to, for example, a pair-factorized state in which there is factorization over pairs 
of adjacent sites in which the stationary probabilities take the form (jTj) where g{m,n) 
is the pairwise weight. The factorized stationary state o is recovered when g{m,n) 
factorizes: g{m,n) = a{m)b{n), in which case the single-site weight f{m) = a{m)b{m). 

Such pair-factorized stationary states have been considered in models |T^ [T6] with 
hopping rates M(mj_i, m*, ruj+i) which depend on the state of both (left and right) 
nearest neighbours: 




g{mi - 1, mj-i) gipruj - 1, m^+i) 

g{mi,mi_i) g{mi,mi+i) 


(33) 


where g{m, n) is the same two-point weight that appears in the expression for the steady- 
state probability (jl]). It has been shown that a pair-factorized stationary state may 





Condensation in models with faetorized and pair-faetorized steady states 


9 



Figure 2. Phase diagram for K{x) ^ and p{m) ~ The spatial extension 

of the condensate ~ where arect = (/? ~ 7 )/(/? ~ 7 + 1 ) for the rectangular and 
cTsmooth = (/3 — 7)/(2/3 — 7 ) for the smooth condensate. Dotted lines correspond to 
constant values of a = 0.05,0.1,..., 0.45. Reproduced from m- 


modify the nature of condensation and allow a condensate spreading over a large, but 
non extensive number of sites. For example, if 

g{m,n) = exp [-J\m - n| + {U/2){5mp + (5n,o)], (34) 

the condensate’s shape is a distorted parabola extending to ~ \/Z sites [15]. References 
|161130] have considered a more general case 

g{m, n) = K{\m — n\) \/p{m)p{n), (35) 

where K{m) and p{m) can be arbitrary, sufficiently-fast decaying functions. It turns 
out that when 

K{x) ~ p{m) ~ (36) 

for a,b,(3,^ > 0 , condensation occurs above a certain critical density of particles if 
7 < 1. The shape of the condensate changes from a single-site one, through a rectangular 
condensate, to a parabolic condensate as fd increases from zero to one, and from one to 
inhnity, see the phase diagram in Fig. O The scaling of the width of the condensate 
with L depends on the parameters / 3,7 in a non-trivial way. These results have been 
recently conhrmed numerically HU, with some small discrepancies attributed to hnite- 
size effects. 

8. Pair-factorized states for discrete time dynamics 

The pair-factorized steady state discussed in the previous section assumes a three-site 
hop rate (l33l) . A very interesting question to ask is whether a two-point u{m,n) such 
as the one in the misanthrope process can also lead to a pair-factorized steady state. 

One might have hoped that when the hop rate u{n,m) does not satisfy the 
conditions for a factorized stationary state, there would still be some choices of rates 
u{m,n) which yield pair-factorized states. However it was shown in |27| that this is not 
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the case. Therefore the misanthrope process either has a factorized stationary state or 
the stationary state has some unknown structure. 

Here we consider a generalisation of the misanthrope process similar to the 
generalisation of the zero range process reviewed in Section [3l That is, we consider 
a Misanthrope-like model with stochastic discrete-time parallel dynamics where all 
particles attempt hops at discrete time steps. Our aim is determine how the structure of 
the stationary state is modified. It turns out that we need to consider pair factorization 
for the steady state probability. 

The dynamics we consider generalise the misanthrope process in two ways. First 
masses move at discrete timesteps therefore there can be simultaneous transport of mass 
at different locations. Second, more than one unit of mass can be transferred from site 
i to i -|- 1. 

More precisely, in each time step, some number /ij of particles depart from site i and 
move to site i + 1 with probability ruj+i). Clearly, we must have 0(/i|m, n) = 0 

for /i > m and for conservation of probability we require Y7^=o '^) — 1- Following 

|19) we will refer to (j){fi\m,n) as the chipping kernel. 

One would like to know whether there is a necessary and sufficient condition on 
4>{fi\m, n) for the steady state P{mi, ..., mi) to factorize over pairs of neighboring sites. 
This means that the following equation must be fulfilled: 

L L 

^ + /ij --F/ii+i-/ii) 

i=l j=l 

X g{nii + yUj_i, rui+i yUj+i - fii) . (37) 


The left hand side of this equation gives the weight (unnormalized probability) of some 
configuration of mass in the system given by {mi}. The right hand side gives the 
sum over the weights of possible configurations at the previous time step multiplied 
by the transition probabilities to the configuration {m*}. The sum over preceding 
configurations is expressed as a sum over the masses, /Xj, transferred from i io i + 1. 
Thus at the previous time step site i had mass + Hi — For a stationary state 
the two sides of (|371) must be equal. 

We now claim that a chipping kernel of the form 

u(p)a(m — fi)h(n -|- /i) 


0 (/i|m, n) = 


g(m, n) 


(38) 


where 


g(m,n) = E u(fi)a(m — fi)b(n + /i). 


(39) 


ii=0 


is sufficient for pair factorization provided that some additional constraint is imposed on 
the functions u{fi), a{m), b(n). Expression fl38|) implies a factorization of the chipping 
kernel into a factor u{fi) which depends on the mass transferred, a factor a(m — p) which 
depends on the mass which remains at the departure site and a factor b(n + p) which 
depends on the resulting mass at the destination site. 
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To demonstrate that the form (155]) indeed leads to pair-factorized stationary state, 
let us insert fl38|) into the stationarity condition Eq. fl37|) . Rearranging indices we obtain 

H Hi+i)) 

* {Mi} * 

= JJ I ^ fii 

i \ 

= Wgijrii+i^mi) . (40) 

jTTTrr require 


In other words, for pair factorization under the form 

/ g{mi+i,nii) = 1 . 


(41) 


The most general solution to this equation assumes the form (for a proof see |27| 
Appendix A) 

h{n) 


g{m, n) = g{n, m) 


h(m) ’ 


(42) 


where h{m) is some function which can be determined by inserting n = 0 into the above 
equation: 

o(0,m) 


h{m) = h(0) 
Finally, from flT2|) and 


^(m,0) ■ 
we obtain 


(43) 


g(m, n)g(0, m)g(n, 0) = g(n, m)g(m, 0)g(0, n) . (44) 

Equation fl44|) . in tandem with the definition fl55]) . is the central result of this section 
and gives a sufficient condition for the stationary state of the generalised misanthrope 
process to take a pair-factorized form. It implies conditions on u{fi), a{m), b{n) which 
we explore in the next section. It remains an open problem whether Eq. fl4T|) is also a 
necessary condition. 

8.1. Single particle hopping 

To simplify the discussion we we now focus on the misanthrope-like case when at most 
one particle can hop from a given site in each timestep. For this purpose we take 
m(p) = + At6fj,^i, and (155]) yields 

g{m, n) = a{m)b{n) -|- Ata{m — l)b{n + 1). (45) 

At this stage At is a parameter but as we shall verify later the limit At —)■ 0 reduces to 
the usual continuous time misanthrope process. We will derive now the relation between 
a{m) and b{m) for arbitrary At > 0. Inserting (155]) in condition (154]) leads, after some 
algebra, to 

m - i{,n) _ ^ m - Hn) _ ( 45 ) 


aim) 


a{n) 
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a(m) = a(m — l)/a(m), (47) 

b{m) = b{ml)/b{m). (48) 

The left and right side of eqnation fj46l) are functions of m and n, respectively. In order 
for (H6l) to be valid for any m, n, both sides of the above equation must be equal to a 
constant k. This gives 

m 


b{m) = 


\d{m) 


1 + Atfe(0)a(m) 


or, equivalently, 


d{m) = k- 


6(0) — b{m) 


(49) 


(60) 


1 + kAtb{0)b{m) 

Therefore, a{m) is determined by b{m) and vice versa, up to two parameters k and At. 
The chipping probability (j){l\m,n) can then be expressed as 

6(0) — 6(m) 


0(l|m, n) = kAtb{n) 


(51) 


1 + kAt(b{m)b{0) + 6(?7,)6(0) — b{m)b{n)) ’ 
and because the hopping probability defines the model, all static and dynamical 
properties are fully specified by giving k, At and one of the two functions d{m) or b{m). 
The pairwise weight function g{m,n) is given by Eq. flT5D . with a{m),b{n) calculated 
recursively: 

1 


a(m = 


n 


ait 


K") =n^(*“ 


(52) 

i=l ^ '' i=l 

where we assumed for convenience that a(0) = 6(0) = 1. This assumption is made 
without loss of generality as it only rescales g{m,n) by a constant factor. 

In the limit At —?■ 0, we obtain factorization of the steady state, since fl45|) tends 
to g{m,n) = a{m)b{n) corresponding to a single-site weight f{m) = a{m)b{m). The 
chipping probability becomes 

(p{l\m,n) = Atd{m)b{n)-i-0{At‘^), (53) 


and in the continuous time limit At —)■ 0 this reduces to a hopping rate 

u{m,n) = a(m)6(n), 

whereas condition (|46|) reduces to 

6(0) — 6(m) 6(0) — 6(n) 

a(m) d(n) ’ 

which is equivalent to the condition (IT^ for the misanthrope process. 


(54) 


( 55 ) 
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In this section we shall outline how the above discrete-time misanthrope process with 
factorized hopping probability exhibits condensation above some critical density of 
particles. We shall learn, perhaps not surprisingly, that the conditions on the hopping 
probability (j){l\m,n) are somewhat similar to those for the hopping rate u{m,n) in the 
continuous time case. 

We begin by defining the grand-canonical partition function 

G{z) = (56) 

{mi} 

where the steady-state weights reads 

L 


i=l 


[a(m,)fe(mi+i) + Ata(mi - l)fe(m,+i + 1)] . 


(67) 


2=1 


We now observe that the expression in square brackets in Eq. fl57|) can be viewed as a 
product of two vectors: 


(a(mj) -|- Ata{mi — 1)) 


b{nii+i + 1) 

and hence the steady-state weight can be rewritten as 

= n + Aio(mi - 1)) ( 


(58) 


2=1 


= Tr 


= Tr 


= Tr 


(a(mj) -F Ata{mi - 1)) 


, 2=1 

■ L 

n 

, 2=1 

■ L 

n 

2=1 


b{mi+i) 
b{mi+i -F 1) 


b{mi) 
b{mi 1) 


{a{mi) + Ataimi — 1 )) 


a{mi)b{mi) Ata{mi — l)b{mi) 

a{mi)b{mi -|- 1 ) Ata{mi — l)b{mi -|- 1 ) 


(59) 


where in the penultimate step we have cyclically permuted the vectors under the trace 
and in the last step evaluated the resulting dyadic product. The grand-canonical 
partition function becomes 

G(z) = Tr p(z)'-] , (60) 

where A{z) is a 2 x 2 matrix: 


kl(^) = 


J2^a{m)b{m)z"^ AtJ2^a{m — l)b{m)z'^ 
0 ‘i'm)b{m + l)z™ At Jfm ~ l)b{m + l)z'' 


(61) 
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The above matrix has two eigenvalues 

A-i-(z) = - (^Aii{z) + A22{z) ± \/{Aii{z) — A22{z)y + 4Ai2{z)A2i{z)'J , (62) 

where Aij is an element of (l 6 T]) . Since A+ is always greater than A_, provided that z > 0, 
\nG{z) = LlnA+ for large system size L, and we obtain that, in the thermodynamic 
limit, the density-fugacity relation is given by 


p{z) = z 


a;(^) 

A+(^) 


(63) 


Since p{z) is an increasing function of 2 ;, condensation will occur if p{z) tends to a finite 
value as 2 ; approaches its maximum allowed value. This maximum allowed values is 
- the radius of convergence of A+(z). From Eqs. (I62|) and (IM]) we see that the radius 
of convergence of A+( 2 ;) is the radius of convergence of any of the four entries of A{z). 
If we focus, say, on Aii{z), we see that condensation criteria reduce to the convergence 
properties of the sum '^^a{m)b{m)z'^. Hence, condensation is possible if 

f{m) = a{m)b{m) (64) 

behaves in one of two ways described earlier i.e. here f{m) plays the same role as the 
single-site weight f{m) in section [T] cases I, II. 


9. Discussion 

In this paper we have given a short review of condensation in factorized and pair- 
factorized states. In particular we have discussed the misanthrope process where the hop 
rates M(mj,mj_|_i) depend on the occupancy of both the departure site i and destination 
site i 1. This process provides a new route to the standard condensation scenario for 
the case of increasing hop rate u{m, n) ~ m'^n^ for 7 > 2 . 

We have studied a generalisation of the misanthrope process to a parallel discrete 
time updating scheme in which simultaneous transfer of mass at different locations can 
occur. We have shown that conditions exist for the stationary state to take a pair- 
factorized form. These conditions generalize the conditions fll2lll3p for factorization 
in the continuous time case. Thus the continuous time factorized stationary state is 
modified into a pair-factorized state when discrete time dynamics is considered. This 
contrasts with the zero-range process in which a factorized form is maintained under 
discrete time updating. 

A major open question which remains is the structure of the misanthrope process 
stationary state for general hopping rates. Also it would be of interest to generalise the 
conditions for factorization and pair factorization in the misanthrope process to more 
general geometries than the one dimensional periodic chain. 
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